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Abstract 

The branching processes in parton showers are studied in perturbative QCD 
for both quark and gluon jets. The emphasis is on the nature of fluctuations 
of both the parton multiplicities and the spatial patterns of the final states. 
Effective measures of such fluctuations are calculated from the data obtained 
*^ I by Monte Carlo simulations. The entropy indices are used to characterize 

Tjlj- ' chaoticity. Both running and fixed couplings are considered. The fixed coupling 

O . case is used to study the onset of chaos. Implications of the results are discussed. 

p^. 1 Introduction 

The possibility of chaotic behavior of particle production in branching processes has 
^ ■ recently been investigated [jl], with emphasis on the search for appropriate measures 

of chaoticity. Unlike classical nonlinear dynamics where the Lyapunov exponent A 
^ i can be defined to characterize divergent distances between trajectories, no such con- 

^ [ ventional description of chaotic behavior of branching processes exist. One of the 

measures found in Q is the exponent k characterizing the exponential dependence 
of the normalized multiplicity variance on the average multiplicity; another is the 
entropy index fig. With those two measures it is possible to distinguish the properties 
of particle production of non-Abelian from Abelian dynamics to the extent that one 
may regard the former as chaotic and the latter not. 

There are two aspects of that investigation that need further extension and ex- 
ploration. The first is an obvious one: the pure gauge theory of gluons only in the 
perturbative QCD branching should be extended to include quarks. The other is 
less obvious. If the QCD dynamics is indeed chaotic, there is a question rooted in 
the conventional chaos theory that should be asked. That is, at what value of the 
control parameter can one identify the onset of chaos? In QCD there is no externally 
controllable parameter with which one can tune the dynamics. Since the nonlinearity 
of the dynamics that gives rise to branching can be represented perturbatively as a 
self-coupling term whose strength is parametrized by the strong coupling constant Og, 
we may regard as the control parameter. Thus in order to investigate the nature 
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of the onset of chaos, we shall reconsider the QCD branching problem by treating as 
as fixed, and then vary it by hand to study the dependences of n and fig on a^. 

In Sec. 2 we treat the full QCD problem with quarks and gluons and with running 
as(g^). In Sec. 3 the onset of chaos is studied by varying the fixed a^. 

2 QCD Parton Showers 

As in [0 0, we follow Odorico's procedure to develop the algorithm for Monte 
Carlo simulation of parton shower in QCD. We shall not repeat here the details of 
the algorithm already described in . What is new now is the inclusion of quarks in 
the branching processes, and the consideration of quark jets in addition to the gluon 
jets. 

The splitting functions are 
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and the running coupling constant is 



"^^^ ' ~ llN,-2Nf ■ MWA^ ■ 

We shall choose massless quark number Nf = 3, color number = 3, and set the 
QCD scale at A = 250 MeV. The approximations made to treat the infrared and 
collinear divergences result in the Sudakov form factors 0, ^ that have the following 
forms for the gluon and quark (or antiquark) vertices, respectively, 

AG(Q^g') =exp|-^fa l'"dzlPo^aa(z)+3Pc^,^(z)]\ , (5) 

A,(Q^/0=exp{-?ft,(||;^)/'"<^.P,_,o(.)} , (6) 

where e = Ql/Q^. We take the end of branching to occur when < Ql = IGeV^. 

Since a gluon can now go into a number of channels, the simulation of a branching 
into resolvable partons requires the designation of a specific final state. To that end 
we calculate the ratio 

^ ^ Ie~' PG^GG{z)dz 

jl-' [Pg-.Gg{z) + ^PG-,qq{z)] dz ' ^ ^ 
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and compare it with a random number G [0, 1]. In < i?, then the G GG 
channel is chosen. Because of the divergence of (|l]) at z = and 1, i? is very close 
to 1; it is ~ 0.8 even for e = 0.3. Thus the self-reproduction of gluons is a far more 
dominant process than qq pair creation. By this fact alone we do not expect the gluon 
jet to change too much by the inclusion of quarks. However, the quark jets are very 
different from the gluon jets, as we shall show. 

We first show in Fig. 1 our result on the distribution of maximum number of 
generations, P{imax), for both the gluon-initiated shower (G jet) and quark-initiated 
shower (Q jet) at two different values of initial virtuality Q. The distributions are 
obtained after 10^ samples are simulated. There are less generations of branching in 
a Q jet than in a G jet because Pq-^qdz) has divergence only aX z = 1 (soft gluon 
limit), while Pg^gg{.z) diverges at both z = and 1. Thus the highly reproductive 
gluons are generally produced with low momenta (and consequently low virtualities) , 
resulting in less multiplicities for the Q jet. For the purpose of studying multiplicity 
fluctuations at different generations of branching, it is meaningful to consider only 
the subset of all showers that have the same maximum generations, imax- We shall 
choose imax to be at the peaks of P{imax)- More specifically, we set imax = 6 and 11 
for Q and G jets, respectively, at Q/Qo = 10^, and imax = H and 22 at Q/Qo = 10^. 

For each fixed i < imax there is a distribution of multiplicities rii at that i. That 
distribution P{ni) is shown in Fig. 2 for the case of Q/Qq = 10^ and for the values 
i = 6 and 11. It is clear that {rii) for G jet is greater than that for Q jet for both i 
values. The fluctuations from {rii) are measured by the normalized variance 

{ni) 

In Fig. 3 are shown Vi vs (n,) for the various cases of jet type and Q/Qq. The most 
striking revelation one sees in that figure is that the general features of the Q and G 
jets do not differ by very much. Focusing on Q/Qq = 10^ we see that their Vi values 
reach about the same peak value even though the corresponding (rii) are different for 
the Q and G jets. That small difference, however, accounts for a difference in the 
value of K defined by 

V, oc (n,)^ . (9) 

Straightline fits to the linear portions for the log-log plots for the Q/Qq = 10^ case 
yield 

0.62, (Q jet, 10 < {ui) < 50) 
~ 0.29, (Gjet, 20 < (ui) < 70) 

This difference in the values of k distinguishes the Q jet from the G jet. It is a 
quantitative measure of the differences exemplified by the thin and thick solid lines 
in Fig. 2. The dotted lines in Fig. 3 represent the results of pure gauge theory with 
Nf = 0. Note that they are indistinguishable from the new results for G jet with 
quarks included. 
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The above result may be regarded as the characteristics of the temporal behaviors 
of the branching processes, when {rii) is interpreted as an analogue of the time elapsed 
in a classical tragectory 0. For a more effective description of the fluctuations, we 
now examine the spatial patterns in phase space at the end of branching and study 
how they fluctuate from event to event. In terms of the final-state momentum fraction 
X of a parton, where x = HiZi, the inclusive distribution p{x) is highly nonuniform 
in the interval < x < 1. Even in terms of the ( variable, where ( = —logi^x, the 
distributions appear Gaussian, as shown in Fig. 4. Apart from the heights of dn/dC, 
it appears that there is little difference between the Q and G jets. However, those 
are only single-particle distributions. To see spatial patterns it is necessary to study 
the normalized factorial moments Fq in small bins. To that end we must first use 
a spatial variable X in terms of which p{X) is uniform. Adopting the cumulative 
variable [||, 0, |, ^ defined by 

^(0= f^PiOdC/ t%{OdC , (11) 

where C1.2 are the extrema of the C, range in Fig. 4, we have p{X) = constant for 
< X < 1. It is this unit interval in X space that we partition into M bins of width 
6 = 1/M. 

For an event e let rij be the multiplicity of partons in bin j at the end of branching 
so that the factorial moment for that event is 

M 

fl{M)=M-'Y.^,{n,~l)...{n,-q + l) . (12) 
The spatial pattern for the eth event at resolution 6 is then described by 

^; = /,7(/iT ■ (13) 

The fluctuation of from event to event is a dynamical property of QCD branching 
that results in the loss of information about the final state. Thus it is important to 
quantify the degree of that fluctuation, which we do by taking the moments of F^: 

where (■ ■ ■) means averaging over all events. The exponent p can be any real number, 
not necessarily an integer. Since F^ deviates more from 1 when 6 is smaller, we search 
for the power-law behavior in M, i. e. , 

Cp,g(M) OC M^^^P) . (15) 

The entropy index is defined by |l], || 

(16) 

p=i 
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The simulated data we use for this part of the analysis are not restricted, as was 
done for the study of Vi where the maximum number of generations, imax, was chosen 
at the peak of P{imax)- In the study of spatial patterns that change from event to 
event, we consider all events whatever imax may be. After simulating 10^ events we 
calculate Cp_g(M), the results for which are shown in Fig. 5 for Q and G jets at 
Q/Qo = 10^. Evidently, those moments are much larger for Q jet than for G jet; 
thus the former has much larger spatial fluctuations. The entropy indices provide a 
more efficient way of describing that property, as shown by the plots in Fig. 6. The 
values of iiq for the G jet are not significantly different from the ones determined in 
[0, where Nf = 0. But they are much smaller than fig for the Q jet. This difference 
is one of the important findings in this work. It could not have been inferred by 
examining the single-particle distributions such as those shown in Fig. 4. Because of 
the differences between Pq^qc and Pq^gg, the evolution, the multiplicity, and the 
distribution of Fq are all different for the Q jet, as compared to the G jet. If the 
branching dynamics for the G jet is chaotic, then it is even more so for the Q jet. 

3 Varying the Fixed Coupling 

As we have found in |^, the non-Abelian QCD dynamics is different from the 
Abelian x model in that k, > and fiq are larger. Whether or not it means that the 
former is indeed chaotic is not certain, since we lack a definitive criterion for chaoticity 
in problems where no trajectory can be defined unambiguously and where the number 
of degrees of freedom increases with evolution. One way of investigating further the 
nature of possible chaotic behavior in QCD is to examine the onset of chaos. Since 
controls the strength of the nonlinearity in the problem, it is reasonable to propose 
that as be varied in order to find the threshold of the chaotic behavior. Of course, it 
means that we must first fix as and calculate the fiuctuation properties of the parton 
showers, and then examine how those properties depend on as- 

With as not running, the Sudakov form factors in and lose one degree of 
log each and become 

AG{Q\K')=exp!^-^en^j'''dz[PG^GG{z) + 3PG^qq{z)]^ , (17) 

Aq{Q\K') = expl^-^in^l'''dzPq_.qGiz)'^ . (18) 

Apart from this change the simulation of parton showers is just as in Sec. 2. We shall 
fix Q/Qq at 10^, and let as vary over the range 0.05 to 0.3, which roughly covers the 
running range as the virtuality degrades from Q to Qq. At the low end of as (near 
0.05) the corresponding value is unrealistic high, but is nevertheless considered 
here because of our interest in the problem at the very small as- 

For fixed as the distributions of imax have shapes that are similar to those in 
Fig. 1 for running as- The peak of P{imax), denoted by imax^ increases with as for 
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both Q and G jets, as shown in Fig. 7. It is important to recognize that the number 
of generations in a branching process is very small at low a^. That is because the 
survivability of a parton without emitting a resolvable gluon or creating a qq pair 
is nearly 1 at very small a^, as is evident in (|l^ and ([18|). That is quite unlike 



the running as{q^) case, where Ois{q^) always gets bigger as q^ evolves toward Q^. 
The smallness of i^ax small as is particularly acute for the Q jet, since when a 
gluon is emitted by a quark, the gluon momentum and virtuality are suppressed, 
thus depriving the Q jet a rich source of parton generations. This property of low 
multiplicity at low will be recalled later to explain other features to be uncovered. 

From the simulated events we calculate Vi as before and find the dependence on 
{ui) as shown in Fig. 8. Evidently, there is a universality in how Vi depends on (nj). 
The increase of merely extends the range. There is a bend in the rise of Vi for G 
jet, but the rise seems to persist for the Q jet, a phenomenon that we have already 
encountered in Fig. 3. From Fig. 8 we gain no insight on the development of chaos, 
if indeed the temporal behavior is chaotic, since there is no value of ag that we can 
identify as the onset of chaos. 

Next, for the fluctuations of spatial patterns we show Cp^q{M) in Figs. 9 and 10 
for Q and G jets, respectively. The general features are similar to Fig. 5 for running 
as{q^)i but there are important differences. It helps to focus on a fixed p and examine 
the dependence on a^. While the scaling behaviors for = 0.15 and 0.25 are rather 
similar, in some cases even hard to distinguish, the behavior for = 0.05 stands out 
markedly different. That is, is much larger at small than at larger as- This 

feature is most pronounced for G jet at q = 2, where the only curves significantly 
different from are for = 0.05. Comparing the vertical scales of the two figures, 
one can also see that the Q jet has larger magnitudes of log Cp^g than the G jet. These 
features are both related to the fact that the parton multiplicity is low for small 
and especially so for the Q jet. Fluctuations are usually large when the averages are 
small. 

The characteristics of Figs. 9 and 10 are represented in capsule form by the entropy 
indices Hg. In Fig. 11 we show fig for various values of a^, while in Fig. 12 the 
dependences on are shown for fixed q. For > 0.15, fig approaches independence 
on as, but as as decreases toward 0, fig increases sharply. This behavior has not been 
anticipated, and is opposite to what is expected if is to play the role of a control 
parameter that can turn on chaotic behavior as it is increased from a small enough 
value. 

As we have noted earlier, the parton multiplicities are low, when as is small, so 
their fluctuations relative to their averages are very large. That is what fig measures. 
Since chaos implies unpredictability of the final state of a parton shower, our results 
indicate that a parton shower is more chaotic at low as than at high as- Thus our 
conclusion here is that as cannot be treated as a control parameter, useful for tuning 
out chaos. Perturbative QCD is such an intricate nonlinear dynamics that it cannot 
be rendered approximately linear by decreasing the strength of the nonlinear term. 
Without the nonlinear coupling there is no evolution, no particle production, and no 
linear dynamics. 
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4 Conclusions 



Our study of the full QCD branching problem that has both quarks and gluons has 
revealed several properties of the parton showers. The gluon jet with quarks included 
is not very different from the case of pure gauge theory [|l], 0] . But the quark jet 
exhibits more fluctuations (relative to the average) than the gluon jet. The measures 
used to quantify the fluctuations are Vi for the temporal behavior and fiq for the spatial 
behavior. They both possess features that are characertistic of chaotic behavior. 

In the fixed coupling case we have found that the results on fluctuations are 
roughly the same as in running problem, provided that the fixed is not set at a 
very low value, outside the range of the running a^. However, in attempting to learn 
about the onset of chaos we have treated as a control parameter and examined 
the parton showers with ag being allowed to approach a very small value. The hope 
was to see how the chaotic behavior would grow as ag is increased. What we have 
found is that the QCD branching processes are more chaotic at very low a^. Unlike 
classical nonlinear dynamics, the nonlinearity in the non-Abelian gauge theory cannot 
be turned off without causing the system to lose the mechanism for temporal evolution 
through branching. A small amount of nonlinearity leads to large fluctuations relative 
to the mean number of partons produced, resulting in more unpredictable final states. 

The measure by which these properties can be described is the entropy index 
fig. We have seen how fig has emerged as a highly effective description of the degree 
of fluctuations of spatial patterns. It is our opinion that the method of analysis 
employed here is not only suitable for the study of parton showers, but also for all 
problems that involve changing spatial patterns, including the conventional problems 
in classical chaos. 
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Figure Captions 

Fig. 1 Distributions of the maximum number of generations in QCD branching pro- 
cesses for quark (Q) and gluon (G) jets. 

Fig. 2 Multiphcity distributions at i = 6 and 11 for quark and gluon jets. 

Fig. 3 Normahzed variance vs average multiphcity as i is increased from 1 to imax- 
The dotted lines are for pure gauge theory without quarks. 

Fig. 4 Single-particle distribution in C = — logx. 

Fig. 5 Log-log plots of Cp^q vs M for various values of p and q. 

Fig. 6 Entropy indices /ig for quark and gluon jets. 

Fig. 7 The dependence of i^^^x on the fixed coupling a^, where i^"^ is the value of imax 
at the peak of the distribution P{imax)- 

Fig. 8 Vi vs (rii) for various values of fixed as- 

Fig. 9 Log-log plots of Cp^q vs M for quark jets only when is set at three possible 
values. 

Fig. 10 Same as Fig. 9 for gluon jets. 
Fig. 11 fig vs q for various a^. 
Fig. 12 //q vs as for various q. 
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